New results on catalyzed BBN with a long-lived negatively-charged massive particle 
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It has been proposed that the apparent discrepancies between the inferred primordial abundances 
of 6 Li and 7 Li and the predictions of big bang nucleosynthesis (BBN) can be resolved by the existence 
of a negatively-charged massive unstable supersymmetric particle (X~) during the BBN epoch. 
Here, we present new BBN calculations with anX" particle utilizing an improved nuclear reaction 
network including captures of nuclei by the particle, nuclear reactions and /3-decays of normal nuclei 
and nuclei bound to the X~ particles (X-nuclei), and new reaction rates derived from recent rigorous 
quantum many-body dynamical calculations. We find that this is still a viable model to explain the 
observed 6 Li and 7 Li abundances. However, contrary to previous results, neutral X-nuclei cannot 
significantly affect the BBN light-element abundances. We also show that with the new rates the 
production of heavier nuclei is suppressed and there is no signature on abundances of nuclei heavier 
than Be in the X~-particle catalyzed BBN model as has been previously proposed. We also consider 
the version of this model whereby the X~ particle decays into the present cold dark matter. We 
analyze the this paradigm in light of the recent constraints on the dark-matter mass deduced from 
the possible detected events in the CDMS-II experiment. We conclude that based upon the inferred 
range for the dark-matter mass, only X~ decay via the weak interaction can achieve the desired 7 Li 
destruction while also reproducing the observed 6 Li abundance. 



The nucleosynthesis of light elements in the big bang 
is a unique probe of new physics which may have oc- 
curred during the first few minutes of cosmic expansion 
in the big bang. Of particular interest in this work is 
the apparent discrepancy between the inferred primordial 
abundances of 6 Li and 7 Li and the predictions of stan- 
dard BBN. A popular model to resolve this discrepancy 
is the existence of an unstable negatively charged super- 
symmetric particle during the nucleosynthesis epoch [lj 
1 1 31 ] - Depending upon their abundance and lifetime 
such particles can catalyze the nuclear reactions leading 
to enhanced 6 Li [l[ and depleted 7 Li [f| || as required 
by observations. Here we present new calculations based 
upon a substantially improved nuclear reaction network 
for this A~-catalyzed BBN. We solve numerically the 
nonequilibrium nuclear and chemical reaction network 
associated to the X~ particle 7] with improved reaction 
rates derived from recent rigorous quantum many-body 
dynamical calculations @. We show that both the 6 Li 
and Li problems can still be solved. However, contrary 
to earlier speculation [icj |. there is no signature in the 
primordial abundances of heavier nuclides produced by 
this mechanism. 

Also in this work we examine the version of this model 
in which the X~ particles decay into the present dark 
matter. In such models the allowed lifetimes and abun- 
dances can be sensitive to the mass of the dark-matter 
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particle. In this regard the recent results of the Cryogenic 
Dark Matter Search experiment (CDMS II) are of inter- 
est. Possible detected events imply an upper limit on the 
elastic scattering spin-independent cross section between 
the weakly interacting massive particle (WIMP) and the 
nucleon [14j. Based upon this, they have identified an 
allowed parameter region of the WIMP mass of 40 GeV 
< wdm < 200 GeV which is consistent with both the 
CDMS II experiment and the DAM A/LIBRA data. We 
discuss the implication of this mass constraint and show 
that the 7 Li problem can still be resolved together with 
the 6 Li abundance, but only if the negatively charged 
particles decay into a lighter dark-matter particle via a 
weak charged boson exchange. 

The primordial lithium abundances can be inferred 
from measurements of absorption line profiles in metal- 
poor stars (MPSs). These stars exhibit roughly constant 
values of the abundance ratio, 7 Li/H, as a function of 
metallicity (TB-dH implying a primordial abundance of 
7 Li/H= (1 - 2) x 10" 10 . The standard BBN model, how- 
ever, predicts a value that is a factor of 2 — 4 higher 
(e.g., 7 Li/H=(5.24i°;^) x 10" 10 [22]) when one uses the 
baryon-to-photon ratio determined from an analysis (23j 
of data from the Wilkinson Microwave Anisotropy Probe 
(WMAP) of the cosmic microwave background (CMB) 
radiation. This discrepancy requires a mechanism to re- 
duce the 7 Li abundance inferred from BBN. The combi- 
nation of atomic and turbulent diffusion [H, [25[ might 
have reduced the 7 Li abundance in stellar atmospheres, 
but this possibility has not yet been established [261 ] . 

An even more intriguing result concerns the 6 Li/ 7 Li 
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isotopic ratios for MPSs. These have been deter- 
mined [18| and a high 6 Li abundance of 6 Li/H^ 6 x 1CT 12 
has been suggested. This is ^1000 times higher than the 
standard BBN prediction. One should be cautious, how- 
ever, in interpreting these results in that convective mo- 
tion in stellar atmospheres could cause systematic asym- 
metries in the observed stellar line profiles and thereby 
mimic the presence of 6 Li [27]. Nevertheless, several 
MPSs, continue to exhibit high 6 Li abundances even af- 
ter carefully correcting for the convection-triggered line 
asymmetries (28j . 

Be and B abundances have also observed in MPSs. 
9 Be [29T - I34} and B [35T - l38j abundances appear to increase 
roughly linearly as the Fe abundance increases. The ab- 
sence of a plateau in the abundances of Be and B at low 
metallicity, however, suggests that these elements are not 
of primordial origin. 

Nonthermal nuclear reactions induced by the decay of 
exotic particles have been studied (39l . |40| as a means 
to provide a cosmological solution to the Li problems. 
Nonthermal reactions triggered by the radiative decay of 
long-lived particles can produce 6 Li nuclides up to a level 
~ 10 times larger than the observed level without causing 
discrepancies in abundances of other light nuclei or the 
CMB energy spectrum [401 ] . 

Another solution to the lithium problems of particu- 
lar interest here is that due to t he p resence of negatively 
charged massive particles X~ [4lM43j | during the BBN 
epoch. They affect the nucleosynthesis in a different 
way fll-fl3|. The X~ particles become electromagneti- 
cally bound to positively charged nuclides with binding 
energies of ~ 0(0.1 — 1) MeV with the largest binding en- 
ergies for heavier nuclei with larger charges. Since these 
binding energies are low, the bound states cannot form 
until late in the BBN epoch. At the low temperatures as- 
sociated with late times, nuclear reactions are no longer 
efficient. Hence, the effect of the X~ particles is rather 
small. Interestingly, however, the X~ particles can cat- 
alyze the preferential production of 6 Li [1| along with the 
weak destruction of 7 Be [!, 0] . 

A large enhancement of the 6 Li abundance was first 
suggested [l[ to result from an A~ bound to 4 He (de- 
noted as 4 Hex)- This enables the A~-catalyzed transfer 
reaction of 4 Hex(d,A~) 6 Li, whose cross section could be 
seven orders of magnitude larger than the corresponding 
BBN 4 He(c?,7) 6 Li reaction. The cross section for this re- 
action, however, was calculated in a more rigorous quan- 
tum three-body model [J] and shown to be about an or- 
der of magnitude smaller than the estimate adopted in 
Ref. Q. 

Additional enhancements in A "-catalyzed transfer re- 
action rates for the 4 He x (i,A~) 7 Li, 4 Re x ( 3 He,A-) 7 Be, 
and 6 Lix (p,A~) 7 Be reactions were assumed in Ref. 
The rates for those reactions are, however, not as greatly 
enhanced as that of the 4 Hex(d,X~) 6 Li because they 
involve a Al = 1 angular momentum transfer and con- 
sequently a large hindrance of the nuclear matrix ele- 
ment Q. This has been confirmed in recent detailed 



quantum many-body calculations @. 

The resonant 7 Bex (p,7) 8 Bx reaction through the first 
atomic excited state of 8 Bx was suggested [3as a means 
to reduce the primordial 7 Li abundance [51] . A rate 
for this reaction has been calculated in a rigorous quan- 
tum three body model [§], which roughly reproduces the 
value of Ref. [5j but is somewhat inefficient in destroy- 
ing 7 Bex- The resonant reaction 7 Bcx+p — ^ 8 B*(1 + , 
0.770 MeV)x — ^ 8 Bx+7 through the atomic ground state 
of 8 B*(1+, 0.770 MeV) x , i.e., an atom consisting of the 
1 + nuclear excited state of 8 B and an X" particle has 
also been proposed | as a process for 7 Bex destruc- 
tion. From a more realistic estimate of binding energies 
between nuclides and X~ particles Q, this resonant re- 
action was found to exist, but the resonance energy level 
is too high to efficiently destroy 7 Bex- 

The s Be x +P -> 9 B3f -s- 9 B x +7 reaction through the 
9 B X 1 atomic excited state of 9 Bx has also been stud- 
ied 0. However this reaction was found to be not 
operative because its resonance energy is relatively 
large (see Table 2 of Ref. 0). A resonant reaction 
8 Bex(n,A~) 9 Be through the atomic ground state of 
9 Be*(l/2+, 1.684 MeV) x , i.e., an atom composed of the 
1/2+ nuclear excited state of 9 Be and an X~ particle, 
has also been suggested as a possible reaction to pro- 
duce mass 9 nuclides [lfj. Kamimura et al. :9|, how- 
ever, adopted a root mean square charge radius for 8 Be 
of 3.39 fm as a more realistic input. They then found 
that 9 Be*(l/2+, 1.684 MeV)x is not a resonance but a 
bound state located below the s Bex+n threshold. The 
resonant s Bex(n,X~) 9 Bex reaction is thus not likely to 
contribute. 

Neutral A-nuclei, i.e., px, dx and tx have also been 
suggested @ as a means to produce and destroy Li 
and Be through two a-induced A~ stripping reactions 
dx(&,X~) 6 Li and tx(a,X~) 7 Li, and three px induced 
stripping reactions px ( 6 Li, 3 Hea) X~, ( 7 Li, 2a) X~ and 
px{ 7 Be, s B)X~ . The result, however, relies on reaction 
rates calculated within the framework of the Born ap- 
proximation, which is a poor approximation in this low- 
energy regime [^, [42J. The rates for those reactions 
and those for charge-exchange reactions of px{a,p)ax, 
dx{ot,d)otx and tx(a,t)ax have been calculated in a 
rigorous dynamical quantum many-body treatment in 
Ref. @. They found that the cross sections for the 
charge- exchange reactions were much larger than those of 
the nuclear reactions, and that the neutral bound states 
Px, dx and tx were immediately changed to ax before 
they could react with ambient nuclei. The late time pro- 
duction and destruction of Li and Be, therefore, do not 
significantly affect the BBN as shown in this letter. 

The solution to the lithium problems in this catalyzed 
BBN model have been explored by solving the full Boltz- 
mann equations for the recombination and the ioniza- 
tion of nuclides and X~ particles coupled to the nuclear 
reactions Constraints on specific supersymmetric 

models through the catalyzed BBN calculation have been 
obtained in Refs. 0, M, 111, H E3 . 
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Candidates for the leptonic X~ particle of interest in 
these models arc the spin supersymmetric partners of 
the standard-model leptons. Such X" particles (and 
their antiparticlcs X + ) would be produced copiously in 
the hot early universe and subsequently annihilate. Their 
annihilations, however, would frcczeout at some epoch. 
The residual X + particles do not affect BBN because 
they do not bind to the positively-charged nuclei. It is 
possible, however, that the decay of both X ± particles 
affect the hnal light element abundances through elec- 
tromagnetic and/or hadronic showers. Here, however, 
we only consider the X-nuclear reactions, and not the 
effect of subsequent decay. 

The binding energies of nuclei bound to X~ particles, 
i.e., X-nuclei, have been derived by taking account of the 
modified Coulomb interaction with the nucleus 0] under 
the assumption that the mass of the X~ particle is much 
heavier than the nucleon mass. 

We performed the detailed network calculation of the 
catalyzed BBN Q taking into account both the recom- 
bination and ionization of X~ particles with nuclei and 
thermonuclear reactions and /3-decay of normal nuclides 
and X-nuclei. We adopt all of the new reaction rates 
from the rigorous quantum many-body dynamical cal- 
culations of Ref. Q. For the 7 Bex(p,7) 8 Bx resonant 
reaction through an atomic excited state of 8 ^x, rates 
for different masses of X~ have been published [9]. For 
our purposes we adopt their rate for an infinite X~ mass. 
Our results are thus completely different from previous 
studies without the use of the new cross sections. 

We adopt the constraint on the primordial 6 Li abun- 
dance from the observations in MPSs. The primordial 
abundances of 6 Li and 7 Li could be higher than the ob- 
served abundances 18] considering the possible effect 
of stellar depletion of initial surface abundances. The 
observed abundances should, therefore, be considered a 
lower limit to the true primordial ones. Since 6 Li is more 
fragile to nuclear burning than 7 Li (24|, its depletion 
factors could be larger than those for 7 Li. We adopt a 
conservative limit of a factor of 10 above the mean value 
of ( 6 Li/H) MPS = (7.1 ± 0.7) x 10" 12 Q3, and a 3 a lower 
limit to the mean value times a factor of 1/3. The limit 
on the 6 Li/H abundance are thus 1.7 x 10~ 12 < 6 Li/H 
<7.1xl0- n . 

Figures QJ, and [TJd show the results of a catalyzed 
BBN calculation. For these figures the X~ abundance 
was taken to be 5% of the total baryon abundance, i.e. 
Yx = n>x/n>b = 0.05, where nx and rib are the number 
densities of the X~ particles and baryons, respectively. 
Figure [1^ shows the evolution of normal nuclei while Fig- 
ure [T|b corresponds to X-nuclei. 

The abundances of the normal nuclei are very simi- 
lar to the standard BBN abundances until the tempera- 
ture reaches Tg ~ 0.5. The X~ particles then combine 
with 7 Be at Tg ~ 0.5 and subsequently 7 Li at Tg ~ 0.3. 
The 7 Bex produced by these X" captures (Fig. [lb) is 
then destroyed by the 7 Bex(p,7) 8 Bx reaction, primarily 
through the atomic excited state of 8 Bx |>j, and sec- 
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FIG. 1: (color online). Calculated abundances of normal nu- 
clei (a) and X-nuclei (b) as a function of Tg (solid lines). The 
abundance and the lifetime of the X~ particle are set to be 
Yx ~ nx /rib = 0.05 and tx ~ oo, respectively. The dashed 
lines correspond to the standard BBN case. 



ondarily through the atomic ground state 8 B*(1 + , 0.770 
MeV) x composed of the 8 B*(1+, 0.770 MeV) nuclear ex- 
cited state and anl" particle @. We have assumed that 
S B X inter-converts to 8 Be*(2+,3 MeV) x by /3-decay with 
a rate given by the normal 8 B /3-decay rate multiplied by 
a correction term (Qx/Q) 5 , where Q and Qx are the Q- 
values of the standard /3-decay and that of /3-decay for X- 
nuclei 17]. The produced 8 Be*(2+,3 MeV)* then imme- 
diately decays to the three-body channel a+a+X~ [?|. 

When the temperature decreases to Tg ~ 0.1, the X~ 
particles bind to 4 He. Then, the ^"-catalyzed trans- 
fer reaction 4 Hex(d,^~) operates to produce normal 6 Li 
and ^\Ax (after the recombination). Because of the small 
binding energies to the X~ (see Table I of Ref. [7|), neu- 
tral X-nuclei do not form until late times corresponding 
to Tg - 0.03 (for t x ), T 9 ~ 0.02 (for d x ) and Tg - 0.01 
(for px)- The neutral X-nuclei then mainly react with 
4 He nuclei to lose their X~ and to produce 4 Hex (as a 
result of the precise calculation Q) so that abundances 
of neutral X-nuclei are kept low. Nuclear reactions trig- 
gered by neutral X-nuclei are thus not important. 
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FIG. 2: (color online). Contours of constant lithium isotopic 
abundances relative to observed values in MPSs, i.e., d( 6 Li) 
= 6 Li Calc / 6 Li° bs (solid curves) and d( 7 U) = 7 Li Calc / 7 Li° bs 
(dashed curves). The adopted observational abundances are 
7 Li/H= (1.23t°H) x 10_1 ° 01 and 6 Li/ H = C 7 - 1 ± °- 7 ) x 
10 -12 [ll|. Thin solid and dashed lines around the lines of 
d( 6,7 Li) = 1 correspond to the 1 a uncertainties in the observa- 
tional constraint. The gray region is observationally excluded 
by the overproduction of 6 Li. 



Figure [2] shows the contours of <i( 6 Li) (solid curves) 
and <i( 7 Li) (dashed curves) in the parameter plane of 
the abundance Yx and the lifetime tx of the X~ parti- 
cles. d( A Li)= A Li Calc / A Li ohs is the ratio of the calculated 
abundance to the observed abundance. The solid curves 
labeled (i( 6 Li)=10 and 0.2 correspond to upper and lower 
limits on the abundance constraint (1.7 x 1CP 12 < 6 Li/H 
< 7.1 x KT 11 ). Above the d( 6 Li)=0.2 line, 6 Li is pro- 
duced at the observed level in MPSs. 

The thick dashed curves are for <i( 7 Li)=l, 2 and 3. 
The dashed curve for e?( 7 Li)=2 intersects the contours 
of d( 6 Li) for Y x > 1 and t x « (1 - 2) x 1CT 3 . In the 
region above this curve the 7 Li abundance is lower than 
7 Li/H« 2.5 x 10~ 10 . The updated parameter region for 
a simultaneous solution to the 6 Li and 7 Li abundances 
in BBN with negatively charged particle is: Yx ^ 1 and 
T X ~ (1 - 2) x 10 3 s. 

For Tx ^ 10 4 s and Yx ^ 0.3, the calculated abun- 
dance of 7 Li increases slightly due to the 4 Bex(t,X~) 7 Li 
and 4 Hex( 3 He,X~) 7 Be reactions which produce some 
amount of 7 Li. However, this parameter region is not 
allowed due to an extreme overproduction of 6 Li. The 
gray region in Fig. [5] is the parameter region excluded by 
the overproduction of 6 Li. 

If the X~ particle decays via the weak interaction, 
7 Bex converts to 7 Li by a weak charged current tran- 
sition from X- to X°, i.e. 7 Be x -> 7 Li+X° @,[l2|. This 
case is shown in Fig. 6 of Ref. [jj. The results from 
that study do not change by the implementation of the 
new cross sections adopted here. The larger destruction 
rate associated with the 7 Bex — ^Li+X decay followed 



by the 7 Li(p,a) 4 He and 7 Li(X~ ,j) 7 Lix(p,2a)X~ reac- 
tions [1, 0] or a further conversion of 7 Lix by the weak 
interaction [T2J shifts the contours of the 7 Li abundance 
toward smaller values of Yx- In this case the parame- 
ter region which solves both the 6 Li and 7 Li problems is 
Y x « 0.04 - 0.2 and t x » (1.4 - 2.6) x 10 3 s. 

We now consider a model in which the present cold 
dark matter (DM) was produced by the decay of X ± 
particles, i.e., Y DM > Yx- The WMAP-CMB con- 
straint on the cosmological density parameter of cold DM 
^CDM = 0.2 then corresponds to wdm^dm < 4.5 GeV. 
The constraints on Yx required to resolve the 6 Li and 
7 Li problems then imply a range for the the dark- 
matter mass todm- In the case where the reactions 
7 Be x +P ^ 8 Bf ^ 8 Bx+i| and 7 Be x + P ^ 8 B*(1+, 
0.770 MeV)x -^ S B X +7 destroy 7 Be x , the DM mass 
is thus constrained to be muM < 4.5 GeV. On the other 
hand, when the 7 Be x — ^ 7 Li+V° reaction [fl [l2| is in- 
cluded, the allowed mass range increases to todm < 
20 - HOGeV. 

Comparing this result to the allowed parameter region 
for the DM mass of 40 GeV < m DM < 200 GeV [3, im- 
plies that only anl" particle which decay via the weak 
interaction can have existed with sufficient abundance to 
reduce 7 Li produced from BBN. On the other hand, if 
the X~ particles do not decay via the weak interaction 
they are excluded. 

In summary, we have re-investigated BBN in the pres- 
ence of negatively-charged massive particles X~ by solv- 
ing the rate equations with an improved nuclear reaction 
network code |7| • We have adopted the newest quantum 
many-body dynamical calculations 0]. With the new 
rates, we find that, contrary to the speculation in pre- 
vious studies, the neutral JT-nuclei, i.e., px, dx and tx, 
do not significantly affect the BBN abundances. Further- 
more, based upon constraints for the mass of DM parti- 
cles from the possible CDMS-II events we conclude that 
only X~ particles which decay into the DM via the weak 
interaction can simultaneously reduce 7 Li to the desired 
level while producing enough 6 Li in the early universe. 

Finally, we note that in this revised catalyzed BBN 
model there is no signature in the abundances of nuclei 
heavier than Be. Thus, if one were to find a primor- 
dial plateau abundance of Be or B, it would require an 
origin other than this catalyzed BBN model. There are 
three physical processes by which enhanced abundances 
of nuclei heavier than Be could have been formed. The 
first is the cosmological cosmic ray nucleosynthesis in- 
duced by supernova explosions during an early epoch of 
structure formation [46]. This model can resolve the 6 Li 
problem and leave possible abundance plateaus of 9 Be 
and B [13, SI] . The second is the BBN model including 
a long-lived strongly interacting massive particle. Signa- 
tures of such particles are possibly left on the primordial 
abundances of Be and B which may be found in future 
astronomical observations of MPSs [49|]. The third is the 
inhomogeneous BBN model which can lead to a high pri- 
mordial abundance of 9 Be [501 ] . 
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